Fractional Quantum Hall Effect in the Second Landau Level: the Importance of 
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Exact diagonalization of a two-dimensional electron gas in a strong magnetic field in the disk 
geometry shows that there exists a filling factor range in the second Landau level where the states 
significantly differ from those in the lowest Landau level. We show that the difference arises because 
the interaction between composite fermions is not negligible in higher Landau levels, as indicated 
by a substantial mixing between composite-fermion Landau-like levels, or A-levels. We find that the 
exact ground state is well reproduced by composite fermion theory with A-level mixing incorporated 
at the lowest level of approximation. Using the same variational approach in the spherical geometry 
we estimate the excitation gap at filling 1/3 in the second Landau level (which corresponds to 7/3 
of experiment). 
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I. INTRODUCTION 

A two-dimensional electron system (2DES) in a high 
magnetic field and low temperature exhibits a rich struc- 
ture of phases as a function of the filling factor v — 
phc/eB (p is the electron density and B is the magnetic 
field) . The integral quantum Hall effect (IQHE)i is char- 
acterized by quantized plateaus in the Hall resistance, 
where p xy = h/ne 2 , and a vanishing longitudinal magne- 
toresistance p xx around v k, n. This phenomenon occurs 
because the 2DES at integer filling factors v = n is an in- 
compressible liquid with a finite excitation gap. The frac- 
tional quantum Hall effect (FQHE) 2 » refers to plateaus 
at p xy — h/fe 2 in the vicinity of v rs /, where / is a 
fraction. Our understanding of the physics of the lowest 
Landau level (LLL) has evolved greatly in the last two 
decades. The FQHE has been explained by the compos- 
ite fermion (CF) theorjiiii as the IQHE of quasiparticles 
called CF's, which are electrons bound to an even number 
(2p) of quantized vortices of the many-body wave func- 
tion. CF's feel a reduced magnetic field B* = B — 2pp<fio, 
where <^o = hc/e is the magnetic flux quantum. Compos- 
ite fermions form Landau-like levels in the reduced mag- 
netic field, which we will call A-levels (AL's). A-levels 
are analogous to Landau levels (LL's) of electrons at B*. 
(The A-levels have been called CF-quasi-Landau levels in 
the literature. To avoid confusion with the Landau levels 
of electrons, we prefer to use A-levels. "Landau level" will 
refer below exclusively to electronic Landau levels. Com- 
posite fermions can fill many A-levels within one Landau 
level.) When CF's fill an integral number n of A-levels, 
an incompressible quantum liquid with a finite excitation 
gap results, producing an IQHE of CF's. Such condi- 
tions are achieved at electron fillings v — n/(2pn± 1), 
which are precisely the prominently observed fractions. 
Further, the microscopic wave functions based on this 
physics provide an excellent account of the actual eigen- 
functions. CF theory with neglected inter-CF interaction 
thus successfully explains most of the observed fractions 
in the LLL. 

This paper is concerned with fractional quantum Hall 
effect in the 2 nd LL. Given that more than 50 fractions 



have been observed in the lowest LL, one might expect 
a large number of fractions in the second LL as well, 
but FQHE is relatively scarce in the second or higher 
LL's. However, with improved experimental conditions 
(higher mobilities and lower temperatures), many frac- 
tions have been observed outside of the LLL. As seen in 
the experiments of Xia et a/£, the observed fractions are 
2/W = 1/3,1/5,2/5,2/3, and 4/5, where i/W is the fill- 
ing factor of the second LL. (The total filling factor is 
v = 2 + z/ 1 ) or v = 3 + v^ x \ with each LL contributing 
two to the filling factor, taking the spin degree of freedom 
into account.) Gervais et al& have also seen evidence for 
FQHE at 1/5 and 4/5 in the third LL. Indeed, the de- 
creasing stability of CF formation has been suggested on 
theoretical grounds^, and charge-density wave phases are 
known to be dominant in n > 2 LL'sm. The competition 
between many nearly degenerate ground states is illus- 
trated by the observation of the so-called re-entrant inte- 
gral quantum Hall effect: the system goes back and forth 
between the IQHE and the FQHE state several times 
(see, for example Eisenstein et al~ and Xia et alJi), with 
the re-entrant IQHE interpreted in terms of a localiza- 
tion of a correlated bubble crystal proposed theoretically 
by Koulakov, Fogler, and Shklovski: 8 . 

Because the observed fractions are consistent with the 
expectation from a trivial generalization of the CF the- 
ory to the second LL, it is natural to attempt an explana- 
tion in terms of composite fermions. However, the micro- 
scopic description of the observed fractions in the 2 nd LL 
has not been as successful as in the lowest LL. As noted 
by HaldaneiS and by d'Ambrumenil and Reynolds^, 
a generalization of the Laughlin wave function for the 
ground state at v = 1/3 to the 2 nd LLiS has rather poor 
overlaps with the exact second LL ground state (from 
0.47 to 0.61 for 4 < N < 9 particles), which is to be 
contrasted with near unity overlaps in the lowest LL. In 
fact, exact diagonalization studiesiSiii on small systems 
have been unable to capture conclusively the incompress- 
ibility of v^- 1 ' = 1/3, because the system is compressible 
for some particle numbers and incompressible for others, 
and for the incompressible states, the gap varies widely 
with the number of particlesii. The observation of many 
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FQHE states in the filling factor range 2 < v < 4 has mo- 
tivated us to seek a better quantitative understanding of 
the FQHE in the second LL. 

Another motivation arises, interestingly, from certain 
new FQHE states within the lowest Landau level. Pan et 
a/pLi have observed several fractions (4/11, 5/13, 6/17, 
4/13, 5/17, 7/11) that fall outside the primary Jain se- 
quences v = n/(2pn± 1). Although these sequences ex- 
haust the possible fractions for a system of noninteracting 
composite fermions, the residual interaction between the 
CF's can generate more fractions, in perfect analogy to 
the appearance of the FQHE for electrons because of the 
Coulomb interaction. The new fractions are interpreted 
in terms of the higher A-level fractional QHE of compos- 
ite fermion o 3 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 . In particular, v = 4/11 is re- 
lated to the FQHE at v* = 4/3 of CF's, i.e., 1/3 of CF's 
in the second CF AL. (In sufficiently high magnetic fields 
the system is fully polarized, and electrons or composite 
fermions can be taken to be effectively "spinless." This 
is equivalent to the limit in which the Zeeman energy is 
much larger than the cyclotron energy. The filling factor 
z/ 1 ) in the second LL corresponds to v = 1 + z/ 1 ) for 
spinless electrons, but to v = 2 + i/W for the real "spin- 
ful" electrons of experiment.) It was shown by Chang 
and JaiiJ& that the actual 4/11 state is extremely well 
described quantitatively by analogy to the 1/3 state in 
the second LL, but not to the 1/3 state in the lowest LL. 
That further underscores the need for a better quantita- 
tive understanding of the FQHE in the second LL. 

Wojs and QuinniS have argued that the difference be- 
tween the LLL and the 2 nd LL physics is due to some 
kind of pairing in the 2 nd LL. They study the occupation 
number of various relative angular momentum m chan- 
nels in the exact ground states for 8 < N < 14 electrons 
in the spherical geometry and find that as the filling fac- 
tor is varied (tuned by the monopole strength Q in the 
spherical geometry, which corresponds to 2Q flux quanta 
penetrating the surface of the sphere), there are peaks 
in the occupation of the m = 5 channel at certain flux 
values that they identify with v = 7/3,5/2,8/3. These 
findings are taken as evidence for pairing of electrons in 
the m = 5 channel. Wojs and Quinn thus assign a quali- 
tatively different physics to FQHE in the second LL than 
that in the lowest LL, as reflected in the fact that, in the 
spherical geometry, their paired state at 1/3 in the sec- 
ond LL occurs at the LL degeneracy Na = 3A — 4 as 
opposed to Nd = 3 A — 2 for the 1/3 state in the low- 
est LL. (The LL degeneracy is related to the monopole 
strength Q by A^ = 2Q + 3 in the second LL.) No explicit 
wave functions have been constructed for the conjectured 
paired states which can be compared to the exact wave 
functions. It is noted that the likely explanation for the 
FQHE at v = 5/2 20 is based on the notion of pairing of 
composite fermions2ii22i2ii2i described by a Pfaffian wave 
function of Moore and flead 2 ^. 

Our approach is different. We attribute the same qual- 
itative physics to the odd-denominator FQHE in the low- 
est and the second LL's (the half filled LL state behaves 



very differently in the two LL's) and argue that the quan- 
titative differences arise because of substantial A-level 
mixing in the second LL. The A-level mixing is a sig- 
nature of the residual interaction between CF's. The 
negligibility of AL mixing for the lowest LL FQHE states 
is taken to imply that the CF's are weakly interacting. 
We will see that composite fermions are more strongly 
interacting in a range of filling factor in the second LL, 
although not so strongly as to destabilize all FQHE com- 
pletely. The pairing at 5/2 is already recognized as a con- 
sequence of a weak attractive interaction between com- 
posite fermions. 

The higher LL FQHE has also been investigated by 
Goerbig, Lederer and Morias-Smith 26 using the Murthy- 
Shankar formulation of composite fermions^ 7 -. They do 
not consider AL mixing in their approach. 

We will restrict electrons to the second LL and neglect 
any LL mixing. (In practice, we map the second LL prob- 
lem into an effective lowest LL problem, and work within 
the LLL.) The AL mixing is treated perturbatively by di- 
agonalizing the full Hamiltonian within a correlated ba- 
sis which includes the "unperturbed" CF ground state as 
well as the "unperturbed" particle-hole pair excitations 
across A-levels. This is called CF diagonalization. The 
resulting ground state incorporates AL mixing, and can 
be improved perturbatively by inclusion of successively 
higher energy CF excitons. With minimal AL mixing 
(keeping the lowest energy CF exciton), CF diagonaliza- 
tion produces explicit wave functions which, when tested 
against exact wave function, are found to be excellent 
approximations to the ground states, thus demonstrat- 
ing that AL mixing captures the physics of the second 
LL FQHE. This will be seen to be the case for the entire 
filling factor range where the second LL behaves differ- 
ently from the lowest LL. In particular, for z/ 1 ' = 1/3, 
the overlap with the exact wave function increases from 
0.71 to 0.95 (for six particles in the disk geometry) upon 
lowest-order AL mixing. 

The AL mixing caused by the residual inter-CF inter- 
actions has been studied through CF diagonalization pre- 
viously in various other contexts, for example the quan- 
tum dot states^ 8 -, the Luttingcr liquid at the edge of the 
FQHES, the stability of the FQH liquid state against 
CF excitons 3 ^, and for obtaining improved variational 
bounds for the FQHE ground state energies in the lowest 
LL 31 . It has quantitative, and sometimes even qualita- 
tive consequences. 

The plan of the rest of the paper is as follows. The 
next section shows, by comparison of the exact ground 
states in the lowest, the second and the third Landau 
levels for six electrons on a disk, that there exists a fill- 
ing factor range where they behave differently. Sec. IIIII 
will outline the method of CF diagonalization, which in- 
corporates the inter-CF interactions through AL mixing. 
Sec. IIVI will include the effect of AL mixing at the lowest 
order to obtain modified variational ground states in the 
second LL, and find high overlaps, between 0.94 and 0.98, 
with the exact ground state. We then proceed to use the 
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same method in the spherical geometry to estimate the 
excitation gap at i/W = 1/3 in Sec. IV! 



II. EXACT DIAGONALIZATION 



Basics 



The Hamiltonian for the 2DES is 
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H K +Hj + guB • S; 
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where is the band mass of the electron and e is the di- 
electric constant of the host semiconductor. We use com- 
plex coordinates z = x ~ iy on the plane and a symmetric 
gauge vector potential A = (— By/2, Bx/2, 0). The dis- 
tance will be measured in units of the magnetic length 
Ib = y^hc/eB and interaction energy in units of e 2 /els- 
The angular momentum J z commutes with both Hk and 
Hi, and the eigenfunctions of J z and Hk may be written 
as 
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where the lowering operators are defined as 



a + P 



-, /? = 5 X + £&, = 2d z . 



(2) 



(3) 



(4) 



In terms of the raising and lowering operators, the kinetic 
energy part of the Hamiltonian is given by 



H K 



a^a + — ) //a.', 



(5) 



with w c = eB/rribC, and the angular momentum operator 
is J z = b^b — a^a. The LL index is n, and the angular 
momentum is I — n. The Slater determinants of states 
given in Eq. © are used as basis vectors when Hj is 
diagonalized in a finite subspace of the Hilbert space. 

We will neglect LL mixing throughout this paper, 
which is a good approximation when the interaction 
strength is small compared to the cyclotron energy, and 
also restrict ourselves to states which are maximally spin 
polarized (i.e. fully spin polarized in the topmost par- 
tially filled Landau level); therefore the last term in 
Eq. (the Zeeman term) can be dropped. The ex- 
act diagonalization will be carried out in each L sector 
separately, where L is the total angular momentum of 



electrons, because the Coulomb interaction does not mix 
states with different L. Fixing the total L confines the 
electrons to a disk (hence the name "disk geometry"). 
The filling factor can be tuned by varying L. The fill- 
ing factor is not really a well-defined quantity for a finite 
system, and finite systems do not necessarily represent 
a thermodynamic system with a well-defined filling fac- 
tor. For certain wave functions it is possible to write a 
relation between L and v for a finite N. Laughlin's wave 
function^ for v — l/(2p+ 1) (where p is integer) satisfies 
the relation 
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2p+ 1 



iV\ 1 
2 T 



(6) 



We thus know what L value corresponds to v = 1/ (2p+l) 
in the disk geometry, and it is natural to assume that the 
filling factor decreases monotonically with L. 

Two-body interactions in the lowest LL are character- 
ized by the Haldane pseudopotentials 10 - 33 V m , which are 
the energies of two electrons in relative angular momen- 
tum m state, \ip m ), 



Vrr, 



(ip m \V(zi ~ z 2 )\ip m ) 



rdrV{r)r 2m e- r2 l A 



(7) 
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where V(r) — l/r is the Coulomb interaction, V(q) is 
its Fourier transform, and L m is a Laguerre polynomial. 
The interaction part of the Hamiltonian can be written 
in terms of the pseudopotentials: 



2 2m+1 m\ 
= ( qdqV{q)L m (q 2 )e- q \ 



Hi 



= ££^ 

i<j m 



(9) 



where projects the wave function of the z,jth parti- 
cles into the state of relative angular momentum m. As 
the spatial part of the fermion wave functions is fully an- 
tisymmetric (full spin polarization is assumed), only the 
odd to channels are filled. 



B. Comparing states in different Landau levels 

When studying the state in nth LL we assume that 
the completely filled 0, . . . , (n— 1) LL's are inert, i.e., we 
diagonalizc in the Hilbert subspace of the nth LL. The 
existence of LL raising operators in the planar geome- 
try lets us map this problem to the lowest LL by using 
effective pseudopotentials VZ: 
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(c,insi-*s,)ic,; 



where \ip% 



= (a^) n \ip m ). One can show that 
V n (q) = (L n (q 2 /2)) 2 V(q), 



(10) 
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where V n (q) in the Fourier transform of the effective in- 
teraction for LL n. Then V£ can be calculated in a closed 
form by Eqs. JHJ and IjlOjl in a straightforward manner. 
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1: Overlap of the LLL ground state I^m) with the 
2 na LL ground state |*S>) (top) and the 3 rd LL ground 

(2) 

state ) (bottom) for N = 6. L denotes the total an- 
gular momentum. The overlap 1(^1° |M?ex ) | is almost zero 
for L = 41,46,52,81, and \{9 ( °J\^J)\ nearly vanishes for 
L = 22,34,41,46,52,56,61,66,71,81. The fillings v = 1/3 
and 1/5 occur at L = 45 and 60, respectively. The lines are a 
guide to the eye. 

Once we map the higher LL problem with a Coulomb 
interaction to a LLL problem with an effective interac- 
tion, we can work within the lowest LL, with different 
interactions simulating different LL's. This is what we 
will do. In the remainder of the article, the phrase "the 
ground state in the nth LL" will really mean "lowest LL 
sibling of the nth LL ground state." Of course, the actual 
nth level ground state can be obtained from the lowest 
LL wave function by promoting it to the nth LL using the 
raising operators. With that understanding, the ground 
states I^^S') of different LL's can be compared by cal- 



culating the overlap (^L |^ 



(n)|,T,(n')\ 



(Of course, the real 



nth LL ground state for n > is orthogonal to the LLL 
ground state.) 

Fig- HI shows how the overlap between the second and 



lowest (upper panel) and the third and lowest (lower 
panel) Landau levels varies with L. It is apparent that 
there is an angular momentum range where the ground 
states in the lowest and the excited LL's are quite dif- 
ferent. As one considers higher LL's, the range of this 
region widens. Similar behavior is found for N = 7 and 
N = 8; see Figs. H and 
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FIG. 2: Overlap of the LLL ground state \^ ( °J) with the 
2 nd LL ground state \^J) (top) and the 3 rd LL ground state 
|*S) (bottom) for N = 7. The overlap \{$ ( °J \$£J) | is near 
zero for L = 40, 46, 52, 70, 95 and quite low (« 0.03) for L = 
61, where L is the total angular momentum. The overlap 
\{^ { °J\^J)\ is rather low over a range of L. The fillings 
v = 1/3 and 1/5 correspond to L = 63 and 105, respectively. 
The lines are a guide to the eye. 

There are angular momenta where the ground states 
in different Landau levels are nearly orthogonal. For 



example, for N = 6 I* 



Wut,(0)\ 



0.00024, 0.005 for 



L = 41,81, respectively, and is within numerical ac- 
curacy for L = 46, 52. This near orthogonality is most 
unexpected. We have found that, in most cases, the fol- 
lowing simple explanation can be given for it. As known 
from quantum dot studies^&M, two-dimensional few elec- 
tron systems in a strong transverse magnetic field have a 
tendency to form crystallites. Crystallite structures with 
N = 6 fall into two classes of symmetry: the (1,5) crys- 
tallite resembles a rotating molecule with one particle 
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FIG. 3: Overlap of the LLL ground state l^ia?) with the 
2 nd LL ground state for TV = 8 particles. The overlap 

K*S|*S>| is almost zero for L = 45,50,57,85,101, where 
L is the total angular momentum. The fillings v = 1/3 and 
1/5 occur at L — 84 and 140, respectively. The lines are a 
guide to the eye. 



at the center and five in a ring; the (0,6) crystallite has 
all particles on a ring. Although the density profile dif- 
ferentiates between the two groups, the pair correlation 
function contains more information about the crystallite 
structure. For N = 6 and L = 41,46, 52, 81 the symme- 
try of the crystallites in the two LL's differ, as suggested 
by the "effective" density profiles p eS (r) in Fig. 0] and 
demonstrated by the pair correlation functions g e (r) in 
Fig. O (The superscript "eff" reminds us that the den- 
sity or the pair correlation function are for the lowest 
LL representation of the nth LL states. The "real" den- 
sity or the pair correlation function in the nth LL can 
be obtained by elevating the wave functions to the nth 
LL with the help of the raising operators.) For L = 41 
and L — 46 the crystallite in the LLL has (1,5) sym- 
metry while the one in the 2 nd LL has (0, 6) symmetry; 
for L = 52 and L = 81 the opposite holds. Although 
the majority of cases of almost zero overlap can be un- 
derstood this way, there are exceptions (N — 7, L = 95 
and N — 8,L = 101). We have not explored this issue 
further. 

Outside a range of L, the higher LL physics is similar 
to that in the lowest LL. It was shown earlienii*^ that 
for v < 1/5 the lowest two LL's should show similar 
behavior. The difference at N = 6, L = 81 seems to 
violate this rule, but the system is so small that this 
exception can be attributed to finite size effects. We have 
found no such exception for N = 7 in the L range that 
we have investigated below v = 1/5 (105 < L < 111); 
for N > 8 the total angular momentum corresponding to 
v = 1/5 could not be reached. 
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FIG. 4: Comparison of the density profiles p cB (r) of the 
exact lowest LL and 2 nd LL ground states for the angular 
momentum (L) values for which the overlap between the two 
states is nearly zero. The reason for the superscript "eff' 
is that the density is calculated for the lowest LL sibling of 
the actual ground state. Note the different behaviors near the 
origin. The two Landau levels prefer crystallites with different 
symmetries at these L values. The distance is measured in 
units of Ib, the magnetic length. 
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FIG. 5: The "effective" pair correlation function g eS (r) of the 
ground state for total angular momenta L = 41 and L = 52 in 
the lowest (n = 0) and the second (n = 1) Landau levels The 
distance r is measured from the center of the disk in units 
of the magnetic length Ib- The reason for the superscript 
"eff" is that the pair correlation function is calculated for 
the lowest LL sibling of the actual higher LL ground state. 
Similar differences were found for L — 46 and L = 81. 



III. CF DIAGONALIZATION 

The idea of composite fermion diagonalization is con- 
ceptually quite straightforward. Instead of diagonalizing 
in the full LLL basis space, as is done in exact diagonal- 
ization studies, we diagonalize the full Hamiltonian in a 
restricted basis of correlated states produced by the CF 
theory. The correlated basis is constructed as follows. 
The CF theory2j4 maps interacting electrons at L into 
non-interacting fermions at L* (taking, for concreteness, 
the disk geometry), 



L* = L-pN(N- 1). 



(11) 



We begin with the Slater determinant states <I>^ , which 
are some low kinetic energy states at L* (chosen as ex- 
plained below), and construct Jain's wave functions as: 



(12) 



where $i is the wave function at v = 1 and Vlll projects 
the state to the lowest LL. In the planar geometry, 



*? p =n(*i- z *) 



2p 



(13) 



j<k 



(The form for <I>i in the spherical geometry is given later.) 
The mapping in Eq. (|12|) amounts to a description of the 
2DES in terms of composite fermions, consisting of the 
bound state of an electron with 2p quantized vortices 
of the 2DES wave function. 2 *>CF denotes a CF with 
flavor 2p. The LL's at L* become AL's of composite 
fermions under the above mapping, and the cyclotron en- 
ergy at L* becomes the effective cyclotron energy of CF's. 
While the diagonalization in the basis defined by 'F^ is 
in principle straightforward, its actual implementation 
requires several technically challenging steps which have 
all been demonstrated previously. The lowest LL projec- 
tion will be evaluated following Jain and Kamilla*. The 
states 'F^ are not necessarily orthogonal; their orthog- 
onalization and the diagonalization of the Hamiltonian 
matrix requires the evaluation of large dimensional inte- 
grals. Mandal and JainSi have demonstrated how this 
can be carried through in a Monte Carlo scheme. We 
note that Monte Carlo is used only to evaluate integrals. 

Next we come to the choice of the basis } at L* . 
Electrons at L* will in general occupy several Landau lev- 
els, which map into AL's of composite fermions through 
the above construction. At the simplest level (referred 
to as the zeroth order approximation), these basis states 
are chosen to be the distinct degenerate kinetic energy 
ground states of non-interacting fermions at L* . The AL 
mixing corresponds to LL mixing of fermions at L* , and 
can be incorporated into theory by considering a larger 
basis: 



H*ff J }, 



M J) }] 



(14) 



where { < &^ J ' ) } represents all rj basis wave functions at L* 
with kinetic energy J (in units of the cyclotron energy) 
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relative to the ground state kinetic energy. The corre- 
sponding CF basis is given by 



(15) 



with * 



n 



A diagonalization in this ba- 



sis incorporates the effect of AL mixing perturbatively. 
Fig. |(| explains the physics pictorially. We will denote by 
X^ JS> and E^ the ground state wave function and the 
ground state energy obtained by CF diagonalization. In 
the spherical geometry, E^ J \L) and x (P) w iU refer to 
the ground state energy and its wave function obtained 
by CF diagonalization in the orbital angular momentum 
L sector. As with all variational wave functions the en- 
ergy E{ J > is a strict upper bound on the exact energy. 
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FIG. 6: The CF basis states (schematically) for N = 6 elec- 
tron system at L = 60, which maps into L* = L — 2N(N — 1) — 
of 4 CF's (composite fermions carrying four quantized vor- 
tices). The y-axis shows the AL index, and the x-axis the 
angular momentum index. The top row shows the state with 
the smallest effective cyclotron energy; there is only one state 
at that energy. The remaining figures show all the linearly in- 
dependent configurations with one additional unit of kinetic 
energy. (There are a total of 10 distinct electron states in the 
J = 1 basis at L* , which produce 9 linearly independent CF 
basis states.) 

Of course, by increasing J, one will eventually obtain 
the exact state (although through an unnecessarily com- 
plicated route). The great advantage of the CF theory is 
that a very good approximation to the exact state is ob- 
tained by diagonalizing in a very small basis, in addition 
to providing an understanding of the physics. For many 
total angular momenta L a unique state is obtained at 



the lowest level (i.e. <I>( ) is unique), which gives a vari- 
ational state with no free parameters. In the spherical 
geometry used in Sec. this happens only when v* is 
an integer. On the disk there are many such states (c.f. 
Table [nj, labeled compact states^. 

The edge effects are significant in the disk geometry 
used in this section, while the spherical geometry used 
in Sec. E] represents the bulk. In the spherical geom- 
etry, excitations are obtained by promoting a CF from 
the highest filled AL to the lowest empty AL, leaving be- 
hind a hole. Thus in ^ij J \ J CF particle-hole pairs exist. 
We will find the Jth order variational state in the sub- 
space spanned by CF states with < J holes in the highest 
filled CF AL of $f(°\ In the disk geometry, once a CF is 
raised to a higher AL, its angular momentum is no longer 
fixed, and one can redistribute the angular momenta of 
the particles in the lowest AL. This additional variational 
freedom lets us investigate edge effects as we\\m£^ 

Once the CF basis is constructed, one evaluates the 
overlaps {$?a \Wp } between the basis vectors and the 



matrix elements of the interaction (i&ci 1 ' \Hi\$!^g ) for 
jijj2 < J in the Jth order using a Metropolis Monte 
Carlo technique. That requires the knowledge of the 
real-space interaction corresponding to the effective pseu- 
dopotentials . The inverse Fourier transform of V n (q) 
in Eq. 1)1 Oil is not defined, as the required integral is di- 
vergent. This, however, causes no problem. Because the 
pseudopotentials completely determine the interaction, 
any real-space interaction V n (r) that gives the same 
through 



yn 



(16) 



will be sufficient for our purposes. Thus we can use some 
convenient prescription for the analytic form of V n (r). 
The goal is that the Monte Carlo evaluation of the ma- 
trix elements of V n (r), required by the diagonalization 
in the correlated CF basis, should converge rapidly. In 
this paper we will use the following form: 



V cS (r) 



1 



M 

E 

i=0 



ar e 



(17) 



The form of the real-space interaction in Eq. I|17|) is based 
on the observation that the long-range behavior of the ef- 
fective interaction should approach the Coulomb interac- 
tion, hence all corrections must be short-range. Keeping 
enough number of terms in the sum will give as accu- 
rate a representation of the interaction as desired. For 
our purposes, it is enough to keep the first seven terms 
(M = 6) in the sum. We calculate the first seven odd m 
pseudopotentials V***, V^ s , . . . , VJjf from Eq. (fT7|) sym- 
bolically by Mathematica, and determine c,; (0 < i < M) 
to satisfy 



V' 

* I 



eff 



V 1 
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TABLE I: Coefficients in Eq. 117H to produce the effective 



interaction for the second LL. 

Coefficient Vaiue 

c -50.36588 

ci 87.38159 

c 2 -56.08439 

c 3 17.76573 

c 4 -2.97162 

c 5 0.25132 

c 6 -0.008435 



The coefficients Cj thus obtained, and used in our calcu- 
lations below, are given in Table [I] and the resulting real 
space interaction is shown in Fig. [3 Alternative prescrip- 
tions are available in the literature. Refs. 7 37 38 use a 



Gaussian functional form for the same purpose: 

M 
i—0 

This form, however, is less convenient for our purposes, 
because it yields a real-space potential that oscillates 
strongly as a function of distance, causing an extremely 
slow convergence in Monte Carlo. 







2 nd L 


_ effective 
Coulomb 


interaction 
interaction 














































































0.5 1 1.5 2 2.5 3 

Distance (l B ) 



FIG. 7: The effective real-space potential V (r) (defined in 
Eq. 1170 which simulates the second LL physics in the lowest 
LL. In spite of the deep dip at r < Ib, where Ib is the magnetic 
length, it reproduces the 2 nd LL pseudopotentials to a very 
high accuracy (details in text). 

We comment on the validity of the truncated interac- 
tion, i.e. Eq. I|17|l with M = 6 and the coefficients in 
Table H] First, the physics of the FQHE is governed by 
the first few relevant pseudopotentials. Since our trun- 
cated interaction reproduces correctly, by construction, 
the first seven odd pseudopotentials, it is expected to be 
quite good. Second, V cS is also guaranteed to reproduce 
the pseudopotentials for very large m, where the inter- 
action approaches 1/r. We have found that the largest 



relative error in V m is 5% at m — 37 up to m < 49. (For 
a given total angular momentum L, the greatest relative 
angular momentum is m max = L — ( Ar ^ 1 ); for our cal- 
culations, which go up to L = 60 for N — 6, the largest 
relevant relative angular momentum is 49.) At such large 
m, the pseudopotential is so small that such error is of 
no consequence. As a final test, we have compared the 
ground state of the truncated interaction with the ground 
state of the Coulomb interaction. In the L < 75 range the 
smallest overlap occurs at L = 71, which is 0.998303. To 
quote a typical number, at L = 45 the overlap is 0.99997. 
The truncated interaction is thus essentially exact. 

The overlap and interaction matrix elements for the 
effective interaction were calculated by the Metropolis 
algorithm. Typically 6.6 x 10 7 to 7.5 x 10 7 Monte Carlo 
steps (MCS) were performed, where a MCS is defined 
as a number of iterations during which each particle co- 
ordinate is expected to be updated once. Averages and 
error bars were calculated from five independent runs. 
When a sufficient accuracy was reached, the standard 
Gram-Schmidt procedure was used to find an orthonor- 
mal basis (often telling us that not all of the vectors in 
the correlated CF basis were independent), and H] was 
diagonalized. 

IV. RESULTS 

We have performed CF diagonalization in the 34 < 
L < 60 range. This encompasses v = 1/3 (which occurs 
at L = 45). The flavors 2 CF's are used for L < 48, 
4 CF's for L > 49. (This choice ensures the smallest CF 
basis.) The results are shown in Table |n] and Fig. [8] 
The zeroth-order variational energies and overlaps are 
given for all 34 < L < 60. The first-order wave function 
is evaluated for all L's except for L — 43, 44; for these 
cases the correlated CF basis is too big (D^ = 83 and 
111, respectively) preventing an evaluation of the matrix 
elements of Hi with sufficient accuracy. 

The principal result of our study, as seen in Table ITT1 
and Fig. [SJ is that the CF theory without AL mixing 
is not satisfactory, but very good agreement is obtained 
once AL mixing is incorporated at the lowest order. Let 
us consider a few specific cases. 

At v = 1/3 (L = 45) the zeroth-order variational wave 
function x is the Laughlin wave function, the energy 
of which (E^p = 3.0796(3)) significantly overestimates 
the exact energy 3.05354. The overlap is 0.712(2), which 
is similar to the values that were found in the spherical 
geometryii. The first-order CF diagonalization yields 
Egl = 3.0621(1) and \(^ ex \ X {1) )\ = 0.9467(8). 

We commented earlier that for certain L values the 
ground states of different LL's have anomalously low 
overlaps. Let us consider L — 41, 46, and 52 for N = 6. 
(We could not investigate L — 81 due to computational 
limitations.) At the zeroth order (i.e., without AL mix- 
ing), the overlaps for these L are K^E'exlx' ')! = 0.614(2), 
0.004(2) and 0.933(7), respectively. For L = 46 the small 
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TABLE II: Comparison of the energies (in units of e 1 jtls) of 



the exact ground state and the states obtained from th and 
1 st order CF diagonalization— . The number in parentheses 
is the error in the last digit. D ex , Dq F and Dq P are the di- 
mensions of the full Hilbert space, of the CF basis in the th 
order calculation, and of the CF basis in the 1 st order calcu- 



lations, respectively. The I s order energy was not calculated 
for L = 43, 44 because D^ F is very large, which makes the 
numerical calculation very time consuming. 



L 




F (0) 
^CF 






D (o) 

U CF 


D (i) 

U CF 


34 


3.57207 


3.5801(1) 


3.5748(7) 


235 


4 


27 


35 


3.48739 


3.5144(1) 


3.49786(2) 


282 


1 


10 


36 


3.44351 


3.47805(2) 


3.4476(2) 


331 


2 


18 


37 


3.39199 


3.42079(1) 


3.3962(1) 


391 


5 


31 


38 


3.37222 


3.3899(3) 


3.3770(5) 


454 


9 


47 


39 


3.31020 


3.3403(1) 


3.3174(3) 


532 


1 


17 


40 


3.25712 


3.27476(2) 


3.2616(1) 


612 


2 


26 


41 


3.22851 


3.2647(3) 


3.2330(5) 


709 


4 


41 


42 


3.18755 


3.2206(3) 


3.1933(9) 


811 


7 


59 


43 


3.15253 


3.1761(1) 




931 


12 


83 


44 


3.12831 


3.1379(4) 




1057 


18 


111 


45 


3.05354 


3.0796(3) 


3.0621(1) 


1206 


1 


28 


46 


3.04126 


3.0805(5) 


3.0484(6) 


1360 


1 


39 


47 


3.01108 


3.0561(3) 


3.0235(5) 


1540 


2 


55 


48 


2.97252 


3.02019(8) 


2.9860(8) 


1729 


3 


74 


49 


2.94205 


2.94937(4) 


2.9443(2) 


1945 


4 


46 


50 


2.87713 


2.88772(2) 


2.8805(2) 


2172 


2 


32 


51 


2.87227 


2.8786(6) 


2.8790(2) 


2432 


1 


19 


52 


2.85393 


2.8667(3) 


2.8597(6) 


2702 


10 


65 


53 


2.82087 


2.8375(3) 


2.827(1) 


3009 


5 


44 


54 


2.77835 


2.79440(4) 


2.7881(5) 


3331 


2 


26 


55 


2.72364 


2.737(2) 


2.7289(2) 


3692 


1 


13 


56 


2.72364 


2.7457(2) 


2.7383(2) 


4070 


5 


39 


57 


2.69540 


2.70221(2) 


2.7018(1) 


4494 


2 


21 


58 


2.68363 


2.7054(7) 


2.6979(7) 


4935 


9 


48 


59 


2.64022 


2.65914(2) 


2.6495(6) 


5427 


3 


25 


60 


2.58918 


2.603(7) 


2.6034(4) 


5942 


1 
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overlap follows because the th order CF basis contains 
only one state, which is very close to the LLL ground 
state. In the other two cases, the th order CF basis 
is larger (with 4 and 10 states); the greater flexibility 
results in better overlaps. With AL mixing, the over- 
laps increase to \{^ ex \x (1) )\ = 0.9814(2), 0.9695(7) and 
0.987(7), respectively. 



ESTIMATION OF THE GAP AT v 



Ci) 



1/3 



As was mentioned above, the existence of an excitation 
gap is crucial for the theoretical explanation of the phe- 
nomena of the FQHE. The best estimate so far is from 
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FIG. 8: Overlap of the exact 2 nd Landau level ground state 
with the wave function obtained from CF diagonalization with 
and without AL mixing (solid circles and stars, respectively). 



the exact diagonalization calculation of Morfi£, where he 
obtained a gap of A 7 / 3 w 0.02e 2 /eiB for N = 10 particles 
at i/ 1 ' = i/3 ; w hich gives a measure of the thermody- 
namic gap. In this section we estimate the excitation gap 
at i/W = 1/3 from CF diagonalization, incorporating the 
residual interactions among CF's. 

For the calculation of excitation gaps we switch to the 
spherical geometry, which is good for examining the bulk 
properties. In this geometry, electrons are constrained 
to move on the surface of a sphere and a radial mag- 
netic field is produced by placing a magnetic monopole 
of strength Q at the center, where 2Q(j>o is the magnetic 
flux through the surface of the sphere (2Q is an integer 
according to Dirac's quantization condition.) 

On the sphere the total orbital angular momentum L 
is a good quantum number. The ground state at L = 
is uniform and spherically symmetric and is considered 
incompressible if it is separated by a finite excitation gap 
from the other states. 

The single particle states Y q i m are called monopole 
harmonics— and are given by 



Ygl m (Sl) 



N qlm {-\) l+m e lq ' l 'u q - m v q+m 
- ? N 



(18) 



l-q 

E 

s=0 



q + l 
I + m — s 



x(v*v) l - q - s (u*u) s , 

where q is the monopole strength, I = q + n is the single 
particle angular momentum where n = 0, 1, . . . is the LL 
index, m = —I, — Z+l, . . . , I is the z-componcnt of angular 
momentum, il = (6, (f>) is the position of the electron on 
the surface of the sphere in the usual coordinates, and 
u = cos(f?/2) exp(— i<fi/2) and v = sin (0/2) exp(i0/2). 
The normalization coefficient is 



N, 



qlrn 



'2Z + 1 (l + m)\(l-m)\ 
~kK {l + q)l(l-q)\ 



(19) 
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The degeneracy of the nth A-level is 21 + 1 and the filling 
factor is defined as 



The distance between particles is taken to be the cord 
distance = 2R\uiVj — ViUj\ where the radius of the 
sphere is R — \[Q in units of magnetic length. 

In the spherical geometry the CF wave function is writ- 
ten in the same form as in Eq. I|12|l 

where $ is now a wave function for TV non-interacting 
electrons at monopole strength q. Composite fermioniza- 
tion (vortex attachment) is again accomplished by multi- 
plication of the Jastrow factor which in this geometry 
is written 

= Y[(ujv k ~ v jUk ) 2 r (21) 
j<k 

Again, $i is the wave function of one filled LL. The pro- 
jection Plll into the lowest LL in the spherical geometry 
is a complicated procedure and interested reader is again 
referred to the literature^. We then arrive at a wave 
function ^ describing N electrons at Q = q + p(N — 1). 

To obtain states corresponding to the filling factors 
(y = n/(2pn + 1)) we create electron states at integral 
filling factors n. In the spherical geometry this amounts 
to setting the monopole strength to 2q — N/n — n, the 
value corresponding to a total degeneracy of N. Af- 
ter attaching 2p vortices to each electron and projecting 
into the LLL we arrive at a wave function describing N 
interacting electrons at filling factor v = n/(2pn + 1) 
called (th is wave function equivalently describes 

non-interacting CFs). As in the disk geometry has 
been shown to be spectacularly accurate when compared 
with exact diagonalization results^. On the sphere 
is a uniform state and exists only at L — 0. With no loss 
of generality we work within the z-component of angular 
momentum L z — subspace of the 2L + 1 degenerate 
states in each L channel. 

The low energy excited states if! a ' are calculated as 

*£> = Plll<S>\ p ^ ] , (22) 

where represents the a states at 2q = N/n — n where 
one electron in the filled LL n has been allowed to occupy 
the LL (n + 1) leaving behind a hole. In general ^ is a 
superposition of Slater determinants. There is one state 
for each total angular momentum L = 1, 2, . . . , (N + n 2 + 
n)/n, so no CF diagonalization is required, though it has 
been shownii that the application of Plll annihilates 
the state at L = 1. 

The calculated excitation spectrum E^\L) for i/ 1 ) = 
1/3 for N — 16 is shown by solid circles in Fig. [§] (We 
are showing a system for N = 16 for illustrative purposes 
only; similar behavior occurs for other system sizes.) It 
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L 

FIG. 9: (Color online) The open diamonds represent the 
energy spectra for the J = 2 basis as a function of total 
orbital angular momentum L, i.e. E^ 2 \L), while the blue 
circles represent the energy spectra using the J = 1 basis, i.e. 
E W (L / 0), at filling factor v m = 1/3 for N = 16 electrons. 
The energies have been corrected for the finite sized deviation 
in the density by multiplication by \J p/ pN = \/2Qv/N. The 
contribution from the uniform positive background has not 
been subtracted. 



has an L = ground state (\J>(°) in this case) separted 
from the other states by a gap. It may be tempting to go 
ahead and calculate the excitation gap, but we know that 
V?' ) is a poor representation of the actual ground state 
in the second LL from the studies of small systems in 
the spherical geometryii and in the disk geometry (this 
work and elsewhere^) and suspect that the excited states 
are also equally bad. Thus, we expand our basis for CF 
diagonalization fEci. I15f) to J = 2. 

The open diamonds in Fig. show the spectra when 
the interaction in Eq. (|17fl is diagonalized using the CF 
basis in Eq. O with J = 2 for N = 16 at u& = 1/3. 
Two points are noteworthy. First, as expected, the low- 
est energy in each L sector has been lowered significantly. 
(Many new higher energy states are also created, but they 
are not of interest in this work.) Second, the gap is signif- 
icantly enhanced. Both these underscore the importance 
of AL mixing for the second LL FQHE. Table ITTll gives 
the overlap between a t L — and for a number 
of system sizes. The relatively small overlaps give fur- 
ther evidence that the CF diagonalization in the J = 2 
basis has significantly improved the ground state wave 
function. 

The overall size of the basis diagonalized for J = 1 con- 
sists of (N + n 2 + n)/n + 1 states across the whole spec- 
trum while the size of the basis diagonalized for J = 2 
is approximately 500 states for N = 16 (9 states in the 
L = channel up to ~ 50 states in the L = 10 chan- 
nel.) That should be compared to the dimension of the 
subspace for the exact state which is ~ 10 10 , making it 
inaccessible to exact diagonalization. The overlap and 
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TABLE III: The overlap between the improved variational 
ground state wave function and the Laughlin wave function 
at u w = 1/3 ((x (2) l* (0) }) for N = 10, 11, 12, 13, and 16. The 
number given in parentheses is the Monte Carlo uncertainty. 



N 



Kx (2) l* (0) )l 2 



10 

11 

12 
13 
16 



0.82(1) 

0.81(1) 

0.82(2) 

0.835(3) 

0.77(4) 



interaction matrix elements needed for CF diagonaliza- 
tion used on average 10 7 Monte Carlo iterations. The 
error reported in each value is the standard deviation of 
the mean calculated between, on average, ten separate 
Monte Carlo configurations. Running in parallel on a 
multi-node computer cluster, approximately 7200 hours 
of computer time were utilized in the spherical geometry 
calculations. To obtain the final energy spectra we im- 
plement the Gram-Schmidt procedure to diagonalize the 
interaction Hamiltonian within the subspace spanned by 
the J — 2 basis as it is generally not orthogonal. 

A technical detail regarding the sampling function is 
noteworthy. We choose this function to be one of the 
basis states that we feel closely represents the true state. 
In principle, in the limit of infinite Monte Carlo itera- 
tions the result calculated approaches the exact result 
regardless of sampling function. In practice, a proper 
choice of sampling function can make the Monte Carlo 
converge faster. Using the effective potential in Eq. I|17|) 
we have found that the Monte Carlo error is on the order 
of 10 times larger than the error when using the pure 
Coulomb potential when \E ,< - ^ is used as the sampling 
function. Consequently, we have calculated the spectra 
using many different sampling functions. Since we are 
only interested in the lowest energy in each L channel, 
this energy was chosen to be the energy with the smallest 
statistical error among the energies calculated with the 
different sampling functions. 

We believe that J = 2 CF diagonalization is quantita- 
tively accurate in the second LL, and proceed to calcu- 
late the excitation gap for FQHE at 2/W = 1/3. Fig.HUI 
shows the energy gap A as a function of wave vector k, 
where Mb = L/R. We have calculated the energy gap 
for N — 10, 11, 12, 13, and 16 represented by the differ- 
ent symbols. The energies for different system sizes have 
roughly collapsed onto a single line indicating that we 
are close to the thermodynamic limit. The Monte Carlo 
error is ~ 10 times what one is used to with gap calcula- 
tions in the lowest LL, making the numbers less certain. 
However, many conclusions can still be reached. 

The shape of A is somewhat different than that of 
the lowest LL 1/3 state which shows one roton mini- 
mum. There seem to be two minima in the spectra, at 
Mb ~ 0.5 and Mb ~ 1.25. The value of roton gap at both 
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FIG. 10: (Color online) The different symbols represent 
the energy gap in units of e 2 /eZ_g for systems sizes N = 
10, 11, 12, 13, and 16 as a function of kls = L/R, where R 
is the radius of the sphere and L is the total orbital angular 
momentum. Each energy has been corrected for finite sized 
deviation in the density as in Fig. [5] 



Mb ~ 0.5 and Mb ~ 1-25 is crudely ~ 0.02e 2 /ds in both 
cases. This is consistent with Morf's estimation 2 ^, fur- 
ther confirming that the CF theory with lowest order AL 
mixing is quite satisfactory for FQHE in the second LL. 
The transport gap is identified with the excitation energy 
at L = N since in the thermodynamic limit this state de- 
scribes an infinitely separated CF particle-hole pair; this 
is estimated to be ~ 0.05e 2 /eZs. The theoretical gaps 
are much larger than the gaps measured experimentally, 
however. In Ref. IZ2I the excitation gaps at 7/3 and 8/3 



were measured as A 



7/3 



0.10 K and A 



8/3 



0.055 K, 



respectively. Using the experimental parameters our cal- 
culated gap for v = 7/3 is A ~ 2 K for the roton and 
A ~ 5 K for the transport gap. The theoretical calcu- 
lations quite generally overestimate the FQHE gap; we 
expect that our estimate would scale down once finite 
thickness, LL mixing, and disorder effects are taken into 
account. In this context, it may be noted that the theo- 
retical gap at v = 5/2 estimated from exact diagonaliza- 
tion studies by Morfi£, A 5 / 2 ~ 0.05e 2 /e^, correspond- 
ing to ~ 5 K, is also large compared to the experimental 
gap^ (~ 0.11 K). 

In principle, our method allows one to obtain the ex- 
citation gap for fillings of the form — n/(2pn + 1). 
However, in practice we find that the Monte Carlo error 
is too large for any reasonable number of iterations, per- 
haps due to the nature of the effective potential (Eq. I17|) 
at short distances and the smallncss of the gap. 



VI. CONCLUSION 

We have shown that the residual interaction between 
CF's accounts for the deviation of 2 nd LL FQHE states 
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and the LLL FQHE states in certain range of filling fac- 
tors. This interaction can be taken into account by allow- 
ing for AL mixing, i.e., letting the ground state hybridize 
with CF particle-hole excitations across AL's, which is 
caused by the inter-CF interactions. The lowest level in- 
corporation of AL mixing produces an excellent account 
of the second LL ground states. We also estimate the 
excitation gap at z/W = 1/3 and find that AL mixing 
strongly renormalizes it upwards. 
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